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Abstract In the context of the dispersive model for 
non-perturbative corrections, we extend the leading renor- 
malon subtraction to NNLO for the thrust distribution 
in e'^e~ annihilation. Within this framework, using a 
NNLL+NNLO perturbative description and including 
bottom quark mass effects to NLO, we analyse data in 
the centre-of-mass energy range -^/s = 14 — 206 GeV 
in view of a simultaneous determination of the strong 
coupling constant and the non-perturbative parameter 
QfQ. The fits are performed by matching the resummed 
and fixed-order predictions both in the R and the log-R 
matching schemes. The final values in the R scheme are 
as{Mz) = 0.113ltEJ™28, ao{2 GeV) 
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1 Introduction 

The thrust (T) distribution [1] in e+e~ annihilation is 
one of the most precisely measured event shape observ- 
able at electron-positron colliders [2-8]. Owing to its 
infrared and collinear safety it can be computed sys- 
tematically in perturbation theory. It spans the range 
between the kinematical configuration of two collimated 
back-to-back jets (dijet limit) and a perfectly spheri- 
cal final state. The deviation from a dijet configura- 
tion is proportional to the strong coupling constant 
and thus allows for a precise determination of its value 
at a given scale. Fixed-order QCD predictions for this 
observable, as well as for other three-jet event shape 
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distributions, are known to next-to-next-to-leading or- 
der (NNLO) accuracy [9,10]. These results were used 
in [11-13] to perform precise fits of the strong cou- 
pling constant. Further fits based on the analysis of 
jet rates have been performed in [14-16]. Fixed-order 
calculations give a reliable description of the experi- 
mental data when three or more final state jets are 
created. On the other hand, they diverge in the dijet 
region (where T — > 1) because of large logarithms of 
the form Q;"ln^"(l ~T). In order to obtain a reliable 
theory prediction for the cross section, these logarith- 
mic terms have to be summed to all orders in the strong 
coupling. The resummed cross section is then matched 
to the fixed-order result so as to obtain a prediction 
holding all over the kinematical range spanned by the 
thrust. The matching procedure is to be defined taking 
care of avoiding double counting with the fixed-order 
cross section [17,18]. Resummed calculations to next- 
to-leading logarithmic (NLL) accuracy are available for 
a wide class of event shape observables [17,19-21], and 
various determinations of as relying on these predic- 
tions matched to NNLO {i.e. NLL-^NNLO) were ob- 
tained in [22-26] . The inclusion of the resummation led 
to an important reduction in the scale uncertainty and 
showed the need to go beyond the next-to-leading log- 
arithmic accuracy. For the thrust cross section, the re- 
summation of large logarithms beyond NLL has been 
performed both following effective theory methods [27] 
and using a traditional approach [28] . A resummed pre- 
diction to this accuracy is also available for the heavy- 
jet mass [29] and the jet-broadenings [20,30]. Moreover, 
precise fits of the strong coupling have been performed 
using the effective theory results [31,32]. 

In this paper we present a determination of as based 
on theoretical NNLL predictions matched to NNLO for 
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the thrust distribution. To describe hadronisation ef- 
fects, we extend the dispersive model [33,34] to match 
NNLL+NNLO and perform a subtraction of the lead- 
ing renormalon at NNLO. Bottom-quark mass effects 
are considered to NLO. The theoretical parton level re- 
sults are reviewed in Section 2. In Section 3 we de- 
scribe the inclusion of bottom-quark mass corrections, 
whereas the non-perturbative corrections and the ex- 
tension of the dispersive model to match the perturba- 
tive NNLL-hNNLO is described in Section 4. Section 5 
is dedicated to the determination of ag, the fit proce- 
dure and the results. Finally, Sections 6 and 7 contain 
the comparison to other recent determinations of 
and our conclusions. 



2 Thrust distribution in perturbation theory 

The thrust variable for a hadronic final state in e+e" 
annihilation is defined as [1] 



A E.|Pi-n| 



(1) 



where pj denotes the three-momentum of particle i, 
with the sum running over all particles. The unit vector 
n is varied to find the thrust direction ht which max- 
imises the expression in parentheses on the right hand 
side. In the present paper we will mostly work with the 
quantity r = 1 — T. 

2.1 Fixed-order and resummed calculations 

The differential thrust distribution in perturbation the- 
ory is known at NNLO [9,10]. At a centre-of-mass en- 
ergy Q and for a renormalisation scale it reads 



1 dcr, ^, _ , ,dA, , _2, .(-IB . , 
- ^ {-r, Q) = as (m) (r) + {n) — (t, x^) 

dC 

+ al{ii) — (T,x^) + 0{ai) , 



where we defined 
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and where a is the total perturbative hadronic cross- 
section for e+e~ — >■ hadrons. The perturbative strong 
coupling constant is defined in the MS scheme unless 
stated otherwise. In Section 4 below we will use a dif- 
ferent scheme for a^. 

In theoretical computations it is customary to normalise 

the distributions to the Born cross section (Tq since, 
for massless quarks, the normalisation cancels all elec- 
troweak coupling factors. However, the experimental 



measurement is of course normalised to the total hadronic 
cross section a. 

The normalised thrust cross-section is defined as 



Rt{t) 



if 
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da {t\Q) 
dr' 



e{T-T')dT', 



(4) 



whose fixed-order expansion has the general form 

Rt (r) = 1 + Ri (r) as{^i^) + R2 (t,/!^) aK^^) 

+ Rs{T,f,')al{„'). (5) 

The fixed-order coefficients Ri, R2, R3 can be obtained 
by integrating the distribution (2) and imposing the 
condition i?T(Tmax,Q) = 1 to all orders, where Tmax 
is the maximal kinematically allowed value (rmax — 
0.4275 at ©(a^)). 

In the two-jet region the fixed-order distribution is 
enhanced by large infrared logarithms which spoil the 
convergence of the perturbative series. The convergence 
can be restored by resumming the logarithms to all or- 
ders in the coupling constant. The matched cross sec- 
tion can in general be written as 



Rt{t) = C{as)E{T, as) + D{t, a,). 



where 



C(a«) = l + ^Cj 



(6) 



(7) 



k=l 

00 n+1 



n=l m=l 

Of 

=Lgi{asL) + g2{asL) + —l3ogs{asL) + ... 

TT 

(8) 

where L = ln(l/r). The function gi encodes all the 
leading logarithms, the function 172 rcsums all next-to- 
leading logarithms and so on. The constant terms Ci 
are required to achieve a full N^+'LL accuracy. D{t, as) 
is a remainder function that vanishes order-by-order in 
perturbation theory in the dijet limit r — ^ 0. It is ob- 
tained by matching the resummed cross-section (first 
term in the r.h.s. of Eq. (6)) to the fixed-order one. The 
matching of the two contributions has to be carried out 
taking care of avoiding double counting of logarithms 
appearing in both of them. To this purpose we follow 
the R and log-R schemes discussed in refs. [17,18,28]. 
The matched cross section and the resulting distribu- 
tion have to vanish at the kinematical endpoint Tmax- 
The resummed contribution does not trivially fulfill this 
property which has to be imposed in Eq. (8) by replac- 
ing the resummed logarithms according to 



p V {xlt)p 



+ 1 



(9) 
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where p is the modificd-logarithm parameter, whose de- 
fault value is set to p = 1. Furthermore in the R scheme 
the single logarithm coefficients G„i have to be modi- 
fied as follows: 











1 - 


( ^ )1 
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Eqs. (9,10) also introduce the resummation scale xl 
which accounts for the ambiguity in defining the re- 
summed logarithms. The resummation functions gila^L) 
as well as the coefficients G„m and also acquire 
a dependence on and need to be modified accord- 
ingly [18,28]. 

In view of matching the NNLL resummed distribu- 
tion to the NNLO fixed order prediction (5) using the 
i?-matching scheme, we need to include the logarithmi- 
cally sublcading terms C2, C3 and G31 in the expansions 
(7) and (8). Eq. (6) then produces a reliable prediction 
of the thrust distribution over a broader kinematical 
range. 

The thrust cross section at NLL accuracy was com- 
puted in [17]. The resummation beyond this accuracy 
was first achieved in [27] using an effective-theory ap- 
proach and revisited in [28] , where the full analytic ex- 
pressions for the 0{al) constant term C2 and the coeffi- 
cient G31 were also obtained. The 0{al) constant term 
C3 is currently unknown, and a numerical estimate is 
given in [28] together with the analytic expressions of 
the functions gi{asL). We report a summary of the re- 
summation formulae in Appendix A. 

To compute the functions gi^a^L), we recast the 
cross-section (4) in the dijet limit as 

Rrir) = ^ ^ ^e-^a^(Q^ a,) + 0(t), (11) 

where {Nq = e''"') 

CTjv(Q',a.) =Hil,asiQ)) 

xP(l,a.(yfg))5(l,a.(^))e-(^). (12) 

The radiator R{No/N) encodes the probability to 
emit a dressed gluon {i.e. with radiative corrections) 
with momentum 



k = up + vp + k±, (13) 

where p and p are the momenta of the back-to-back 
quark and antiquark, respectively. Its expression reads 




The fimction A{as) has a soft origin, while B{as) ac- 
counts for the hard coUincar contributions. They both 
admit a perturbative expansion in Us 

^(a,)=^4^«, fi(a,) = ^4^W, (15) 

I I 

and their expressions are reported in Appendix A. In 

the previous formulae we systematically neglected terms 
of order 0(r) which are absorbed in the remainder func- 
tion D{t, as). The integration countour in Eq. (11) runs 
parallel to the imaginary axis to the right of all singu- 
larities of the integrand. For a description of the terms 
in Eq. (12) (each of which admits an expansion in ag) 
we refer to [27,28]. 

From Eq. (14) we see that the integral over u in the 
exponent is regularised by the lower bound Such 
a bound acts as an infrared regulator which prevents 
the strong coupling constant from being evaluated at 
non-perturbative scales (< Aqcd)- This can be seen 
from Eq. (11) where the contour should be set away 
from all the singularities (in particular from the Lan- 
dau pole). Nevertheless, for resummation purposes we 
can set the contour to the left of the Landau singu- 
larity since it would contribute with a non-logarithmic 
effect suppressed by some negative power of the centre- 
of-mass energy scale. Such contributions, together with 
hadronisation corrections, are numerically sizeable at 
low energies and they are relevant to fit the strong cou- 
pling- 
Several phenomenological hadronisation models can be 
found in the literature: analytical [34 37], empirical [38 - 
40] and statistical [41]. Each of them is based on differ- 
ent theoretical assumptions. In our analysis we restrict 
ourselves to the analytical dispersive model developed 
in [33,34]. 

3 Finite bottom-quark mass corrections 

An additional effect to be considered concerns the as- 
sumption of vanishing quark masses we made in com- 
puting the cross section. Such an assumption is not fully 
justified for low energy data samples {e.g. LEPl, PE- 
TRA energies) for which finite bottom quark mass ef- 
fects are relevant at the percent level [42,43]. We in- 
clude bottom mass corrections directly at the level of 
the matched distribution by subtracting the fraction of 
massless 6-quark events {rb{Q)) and adding back the 
corresponding massive contributions. This results in 

- ^(t,Q) = - ( (1 -rf,(Q))^(r,Q)]i„assless + 

a dr a \ dr 

+r6((5)^(r, 0)1 massive ) . (16) 
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Here ^(r, Q)] massless is the NNLL+NNLO matched 

distribution and ^(r, (5)| massive is the NLO massive 
distribution obtained with the parton level Monte Carlo 
Zbb4 [43]. Since the NNLO correction to the massive 
distribution is currently unknown, we replace it with 
the corresponding massless result 



-Ml • 

V' / I massive 



V JImassless* 



(17) 



The 6-quark mass corrections are generated for the con- 
sidered energies. Futhcrmorc, we perform a bin-by-bin 
interpolation of the energy dependence in order to com- 
pute them for any value of Q included in the considered 
energy range. 

4 Non-perturbative corrections 

In the dispersive model [33,34] non-perturbative cor- 
rections arc accoimtcd for by means of an effective cou- 
pling aeff(fc^) which is supposed to be finite in the in- 
frared region down to fc^ — >■ 0. To introduce the non- 
perturbative correction, we use the paramctrisation (13) 
for the momentum of a soft gluon, or an object with 
gluon quantum numbers {i.e. the offspring of a gluon 
decay). The two Sudakov parameters u and v are re- 
lated by t; = {k\ + rv?)lu, where rv? is the invariant 
mass of the parent gluon. The strong coupling is de- 
fined in a physical scheme as the anomalous dimension 
7^ controlling the evolution of the inclusive soft gluon 
emission probability 



dw = 



2Cf du dk'i 
IT u k^^ 



(18) 



The physical expression of the coupling then becomes 

dacsim'^) 
dm^ 
(19) 



f k 
a,{kl) = r(a,) = a,(0) + / dm \^ ^ 

Jo "-J. "T '''' 



where the effective coupling aefi is defined as the log- 
arithmic derivative of the spectral density for the dis- 
persive relation (19). Integrating (19) by parts (notice 
that q:s(0) = acff (0)) we obtain a dispersive relation for 
the physical coupling as 



oisik'i) 



f 

Jo 



dm 



kl 



{m? + klf 



aeff(m^)- 



(20) 



This relation ensures that the strong coupling constant 
as{k\) remains finite at very low momenta, where the 
perturbative picture fails. It can be recast symbolically 
as the sum of a perturbative and a non-perturbative 
components 



a«(fci) = arHfci) + ar(fci) 



(21) 



None of the two terms in the r.h.s. of (21) is separately 
well defined in the infrared region, but their sum is finite 
down to A;^ = 0. 

The perturbative component in Eq. (21) of the phys- 
ical coupling is defined as the perturbative contribution 
to the physical scheme [44] used above to define the 
strong coupling otg. The inclusive-probability picture 
used in this scheme requires that we integrate inclu- 
sively over the invariant mass rn? of the parent gluon. 
Such an approximation is accurate up to NLL and to 
this accuracy the relation between the physical and the 
MS coupling can be extracted from the soft contribu- 
tion to the Sudakov radiator (14) and reads [44] 



PT 

a. = a. 



1 + 



a, 



A(2) 



(22) 



At NNLL, the non-inclusiveness of the thrust plays an 
important role, and it must be taken into account. We 
need to account for it both at the perturbative as well 
as at the non-perturbative level. At the perturbative 
level we can extend Eq. (22) by including the NNLL 
soft contribution in Laplace space (where the thrust 
factorises) A^^^ , which correctly contains non-inclusive 
effects 



a 



PT 



a, 



^(3) 



+ ... 



(23) 



It should be noted that the NNLL perturbative coeffi- 
cient 

^(3) 

might be different for different observables. 
This would lead to a breakdown of the universality of 
the physical scheme for hadronisation corrections at 
this order. 

According to Eq. (21), the effective coupling ofoff can be 
split into a perturbative term a^g and a non-perturbative 
contribution a^^ which is supposed to be a rapidly 
falling function in the ultra-violet region, where the 
running of the physical coupling (20) matches the solu- 
tion of the renormalisation group equation for the per- 
turbative component oi^"^ . The former is responsible 
for the perturbative strong coupling a^^ in Eq. (21), 
while the latter gives rise to the non-perturbative com- 
ponent a^^. Imposing that for large momenta a^^ does 
not produce any power correction leads to the rela- 
tion [33,45] 



f 

Jo 



dmr 



m^XffK) = 0, n=l,2,... (24) 



By means of Eq. (20) we can find the following useful 
relations between the perturbative and non-perturbative 
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components of a^s and the respective coupUngs 



dkl 







a. 



-a, 







/ 

Jo 



dk^a^^'ikl) 



°° dm^ 



-ma: 



NP 

eflF 



(25) 



(26) 



To obtain Eq. (26) we integrated over k\ and expanded 

the result in a power scries of m? . This is motivated by 
the fact that a^ff (m^) is concentrated at small scales. 
Using Eq. (24) it is straightforward to see that only the 
first term of the scries survives yielding Eq. (26). 
In order to match the non-pcrt\irbativc and the pcr- 
turbative couplings we introduce an infrared matching 
scale /i/ > ^QCD such that a^^{k^) is negligible for 
fc^ > /ij, and the coupling is well approximated by 
the perturbative part a^^ . Since we are only interested 
in the leading power correction (~ 1/Q), we limit our 
analysis to the region where r ^ /ij/Q and wc assume 
that the effective coupling is small enough to neglect 
terms of order 0{alg). It is straightforward to see that 
all the terms but R{No/N) in Eq. (12) are evaluated at 
perturbative scales {i.e. the scales at which the strong 
coupling is evaluated are always larger than the match- 
ing scale Hi). This implies that the only contribution 
to the leading power correction arises from the radiator 
R{No/N). In Eq. (14) the latter is expressed in terms 
of the perturbative coupling in the MS scheme and we 
need to redefine it using the full physical coupling (20). 
We first rewrite Eq. (14) using the identity [46] 



-Nu 



- 1 = -0 U- 



iVo 
N 



+ r2{dlnN)dlNO{u 



N 



O 



No 
N 



(27) 



where 

A(€) = ^(l-e-~T(l-j)) 

and the derivative is meant to act on the whole integral 

whose boundaries arc set by the 0{u — Nq/N) function. 
Eq. (14) then becomes 



(28) 



R 



No 
N 



2/ ^(e-^"-l) / ^Aia.iki)) 

Jo " J«2Q2 Kj_ 



-2 

Jo « 



e-^"-l)S(a,(uQ^)) 



+ dinNBia^i^Q^))^ + O « ln"-2 N) , (29) 



where we used the Leibniz's rule to perform the deriva- 
tives of the integrals. It is straightforward to check that 
the term proportional to C2 contributes to NNLL accu- 
racy. However, we see that it is always evaluated at 
perturbative scales, so it does not contribute to the 
leading power correction. We now observe that the in- 
tegral involving the collinear B{as{uQ'^)) function in 
Eq. (29) gives rise to a subleading power-suppressed 
term which scales as l/Q^, so we can neglect it in the 
non-perturbative analysis. The only contribution to the 
leading power correction stems from the double integral 
in Eq. (29), involving the soft contribution A{as{k'^j). 
Making use of Eq. (23) we can recast this contribution 
using the physical coupling a^'^ (to NNLL accuracy) as 

(No \ 2Cf [Uu ufdkl^PT.,2 

X 0{kl - u'^Q^)0{uQ'^ -kl) + ..., (30) 

where the ellipsis stand for the remaining perturbative 
terms in the radiator. Eq. (30) is nothing but the single 
soft emission contribution to the cross section, with the 
coupling being defined in the physical scheme. 
To correctly account for non-perturbative corrections, 
we replace a^"^ in Eq. (30) with the full coupling (20). 
Moreover, since the physical definition of the coupling 
deals with massive gluons, the massless phase space in 
Eq. (30) is modified to take into account the gluon mass. 
This amounts to performing the replacement k"^ 
kj_ + w? in the two ©-functions in Eq. (30), obtain- 
ing 



+ w? dm? 

m') + ... (31) 



kl 



We integrate by parts over m^ in Eq. (31) getting rid 
of the boundary term (3;s(0), and we integrate over k^: 



2Cf r du 



-Nu 



1) 



dm 



aeff(m^)0(m^ - u^Q^)e{uQ^ -m^) + ... 



(32) 

We now split the effective coupling into its perturbative 
(ct^ff) and non-pcrturbative (a^) components. 

For the perturbative term we use Eq. (25) to repro- 
duce the perturbative soft piece of the radiator (30). 
Replacing the pertiirbative coupling a^"^ with its ex- 
pression (23) we reproduce Eq. (29). We then can write 
the radiator as 
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No 
N 



N 



N 



BF'^{jf-,as) is the full perturbative radiator (29) that 



we recast as 



N 



= - In TV/ii (as In TV) - h2 («« In N) 

- ^/3o/i3(aslniV) + 0«ln"-2iV), 

TT 

(34) 



where the functions hi («« In N) are reported in Ap- 
pendix A. The second term in the r.h.s. of Eq. (33) 
R^^{^,as) is the non-perturbative component, given 

by 

Jo 

(35) 

The non-perturbative term requires some more atten- 
tion. We first recall that wc arc working in the approx- 
imation m/Q <C r, which allows one to expand the 
exponential function as e~^" ~ 1 — Nu + neglect- 
ing subleading terms since they give rise to 0{1/Q-^) 
corrections. It is then straightforward to perform the 
integral over u, obtaining 



2Cp f°° dm'- 9x / 



TT 







(36) 



The term proportional to in the round brackets 
leads to a vanishing contribution because of Eq. (24). 
Making use of Eq. (26) we end up with 



R' 



NP / -^0 



N8Cf 



(37) 



To evaluate Eq. (37) wc replace = oig — ct^"^ ■ Intro- 
ducing the mean value of the physical coupling below 



1 f^' 
M/ Jo 



(38) 



and expanding perturbatively ^ to perform the inte- 
gral, we obtain 



(39) 



Using the expression (23) for a^'^ , the quantity At 

amounts to 



III 8Cf 



Q 7r2 

-a^(^^)- (21n^ + 



TT 
92 



AWpo 

«f(/^?.)4(41n^^+4fln^ + l 

+ 



/3i A(^) 
^'^+2^g+Aa)/3o 



(40) 



Note that the NNLO contribution A^^) in Eq. (40) 
is found to be different from what obtained in [12] 
and [23], where two different assumptions for this new 
term were made. The final expression for A^^^^ [28] is 
reported in Appendix A and it consists of the sum 
of two contributions: the observable- independent three- 
loop cusp anomalous dimension computed in [47] and 
an observable-dependent term proportional to the two- 
loop soft anomalous dimension obtained from the O^a^) 
soft contribution to the thrust cross section [27,28,48, 
49]. The latter turns out to give the leading numerical 
contribution to A^^K 

Using Eq. (40) in Eq. (11) results in a shift of the cross 
section by an amount At. This shift encodes the leading 
non-perturbative correction to the thrust cross section. 
The result in Eq. (40) is not yet complete. So far we 
have considered the dispersive model in its inclusive 
form. The non-perturbative effect of the thrust's non- 
inclusivcness can be accounted for using perturbation 
theory by computing the correction to the form re- 
ported in Eq. (37). Since the physical coupling is de- 
fined as the soft emission probability, one can compute 
the corrections due to incomplete cancellations between 
real and virtual contributions as well as to scenarios in 
which the progeny of the massive gluon goes into oppo- 
site hemispheres. These corrections were computed up 
to 0{al) in [50] where it was shown that they amount 
to a multiplicative {Milan) factor A4 

M = l + (1. 575 Ca -0. 104 tif) = 1.490, np = 3, 

(41) 

where we set the number of active flavours to 3 since it is 
only sensitive to low energy soft radiation. In particular, 
the Up factor in (41) is due to a soft gluon splitting into 
a qq pair of light quarks. It is difficult to estimate the 
uncertainty on At. It can not be excluded that higher 
order 0{al) corrections to the Milan factor could be 
as large as 20% [51]. Wc use this value as uncertainty 
on ^A in our analysis. Since the matched distribution 
is given in terms of a binned histogram, the shift (40) 
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Exp. I Q (GeV) I Fit range | N. Pts. | Ref. 
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cannot be straightforwardly performed since in general 
At is not a multiple of the bin width. We then inter- 
polate between different bins with a cubic spline with 
the; "not-a-knot'' condition in order to evaluate the dis- 
tribution at intermediate values of r. The shift (40) is 
then performed directly on the resulting interpolating 
spline. 

4.1 Hadron masses and decay effects 

In deriving the dispersive contribution to the leading 
power correction (~ 1/Q), all offspring particles pro- 
duced in the m,assive gluon decay were assumed to be 
massless. Taking their finite masses into account would 
change the expression of the thrust in the soft region, 
leading to an additional kind of power corrections [52] 
whose leading scaling is In'^ Q/Q, where A is a constant. 
Such corrections are present in experimental data, where 
massive hadrons contribute to the thrust. Moreover, 
a similar power-suppressed term arises from the final 
state momenta reshuffle due to unstable hadron de- 
cays into lighter particles. Unlike the dispersive con- 
tribution, these power correction cannot be expressed 
as the product of a universal non-perturbative quan- 
tity (e.g. ao) and an observable dependent factor, thus 
they break imivcrsality. A direct consequence is that 
they have a completely different impact on different 
event shape observables. An initial detailed assessment 
of hadron mass effects on event shapes has been pre- 
sented in [52], where the constant A was evaluated un- 
der the hypothesis of local parton-hadron duality, and 
assuming massless partons in the perturbative calcu- 
lation, finding A = 4Ca//3o- Bottom quark mass ef- 
fects in the perturbative contribution could potentially 
modiiy this value. A more recent systematic analysis 
of hadron mass effects on power corrections has been 
performed in Ref. [53]. In Refs. [52,53] it is shown that 
the universality of mass-dependent power corrections 
to two-jet event-shape cross sections can be rescued by 
properly redefining the measurement scheme. This al- 
lows one to study systematically the impact of non- 
perturbative mass corrections for different event-shape 
observables. In particular, they were found to lead to 
sizeable corrections for jet-masses and related quanti- 
ties, while their impact on the thrust is very much lim- 
ited [52,53]. Consequently, we do not include hadron- 
mass effects in our study here. 

5 Determination of ots and ao 

In order to measure the strong coupling constant as and 
the non-perturbative parameter ao we fit the QCD pre- 



Table 1 Data set considered for the simultaneous At of 

as and ao- 



dictions to experimental data reported in Table 1, in the 
centre-of-mass energy range 14 GeV < Q < 206 GeV. 
For consistency we consider data sets which have been 
corrected for QED radiation effects. The other data 
sets would require the inclusion of electroweak correc- 
tions [54] into the theoretical description. Therefore we 
limit our analysis to data measured by the ALEPH [2] , 
the L3 [4] and the TASSO [8] collaborations. 

The theory distributions are computed in form of 
binned histograms. We rebin the theory results in order 
to match the data binning, which is different for each 
experiment. We set the upper limit of the fit range to 
T < 1/3 (the leading order kinematical upper limit) ex- 
cluding the multi-particle region where the theoretical 
prediction fails because of the small number of final- 
state partons considered. The lower limit is set two 
bins away from r = /xj/Q excluding the peak region 
where our treatment of power corrections is not accu- 
rate. Fig. 1 shows results for different fits of the strong 
coupling, obtained by shifting the lower bin around the 
default configuration (labeled with 0). For values of the 
lower bin shift below —2 the peak region is included 
in the fit range, which leads to higher values of Us- 
This region is excluded since the assumption t ^ l^i/Q 
made in computing the leading power correction breaks 
down and subleading contributions become numerically 
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0.5 
0.45 
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R scheme I 








mt, = 4.5GeVl 






E 


J I 1 












-2 2 

lower bin shift 



Fig. 1 Check of the fit stability, when the lower bound of 
the fit range is shifted. The quoted error bars represent the 
statistical uncertainty on the single fit. The plots show that 
the chosen default range is in the middle of a stable plateau. 
Shifting the lower bound of the fit to include more bins in the 
far infrared region leads to large deviations from the stable 
value of «». 



important. The fit is performed by minimising the 
function defined as 




1 da 



cxp 



cxp 



1 da 
a dr 



1 da 
a dr 



(42) 



where Vij are the covariances of the distribution be- 
tween the bins Ti and tj . The general form of the covari- 



ance matrix is Vij = Sij 



Eij, where Sij 



"'stat, i^ij 



is the diagonal matrix of the (uncorrelated) statisti- 
cal errors, while Eij contains the experimental system- 
atic convariances. The diagonal entries of En = a^^^^ ^ 
are given by the systematic uncertainty on the i-th 
bin, while we need to make a plausible assumption on 
the form of the off-diagonal elements. We consider the 
minimal-overlap model, which defines Eij as 



Ei. 



(43) 



The minimization is carried out with the TMinuit 
routine distributed with ROOT and the whole analysis 
was implemented in a C++ code. 

To estimate the fit uncertainties, different approaches 
have been proposed in the literature [18,23]. We decide 
to perform a scan of the theoretical parameters involved 
in the calculation as follows: 

• the renormalisation scale fraction (default value 
— 1) is randomly varied between 1/2 and 2, 



• the resummation scale fraction x^ (default value 
XL = 1) is randomly varied between 2/3 and 3/2, 

• the Milan factor Ai is randomly varied within its 
theoretical uncertainty (20%), 

• the default value of the modified-logarithm param- 
eter p = 1 is replaced hy p — 2, 

• the default 5-quark mass (m^ = 4.5 GeV) is replaced 
by TOfc = 4.0 GeV and mf, = 5.0 GeV. 

The random scans for a;^, xl and Ai are performed 
with high statistics for each combination of mf, and p 
( 3 X 10^ fits in total). The final error on as and is ob- 
tained as the envelope of such variations. Moreover, we 
consider two different matching schemes (R and log-R) 
and take the difference between the two central values 
as systematic uncertainty. The quoted central values 
correspond to the fit results with the default parameter 
setting. 

The scatter plots are shown in Fig. 2. Using the 
R scheme leads to a lower central value and a smaller 
uncertainty on as than in the log-R scheme. We report 
the results separately for the R and log-R matching 
schemes: 



scheme 


as{Mz) 


ao(2 GeV) 


xVd.o.f. 


R 
log-R 


n 11Q1 +0.0028 

U.il01_Q QQ22 

n 1 1 Q7+0.0034 
U.ilO ( _o.0027 


t^S8+°-^°^ 

U.OOO_g 

U.OZ+_Q 


1.54 
1.60 



The difference between the central values is con- 
tained into the scan error bands, so it does not affect 
the final uncertainty. The final uncertainties on as and 
ao can be considered to a good approximation as uncor- 
related. This feature is visible from the scatter plot in 
Fig. 2. To check the dependence of the fit upon the 
infrared matching scale /i/ we perform a fit replacing 
the default value /i/ = 2 GeV by /z/ = 3 GeV. The cor- 
responding results for the mean effective coupling and 
for the perturbative strong coupling are: 



scheme 


as[Mz) 


ao(3 GeV) 


xVd.o.f. 


R 

log-R 


n 11 ■30+00028 
U.iiOU_o.o021 

1136+0 0034 

U.ilOU_Q QQ27 


4'?n+0 068 

U.+dU_o.032 

422+° °'5° 

U.'iZZ_g Q2g 


1.56 
1.62 



The as value is not affected by the infrared match- 
ing scale, whilst ao decreases for larger values of ^7. We 
also observe that the relative uncertainty on the non- 
perturbative parameter ao is smaller for fij — 3 GeV. 
This can be explained observing that the uncertainty 
on ao is mostly controlled by the variation of the Mi- 
lan factor A4 . Given the multiplicative relationship be- 
tween the two quantities, a symmetric variation of A4 
does not correspond to a symmetric uncertainty on ao- 
Moreover, it implies that for smaller central values of 
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[25] (NLL+NNLO) 
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This fit 
(NNLL+NNLO) 



0.105 



World averager 
0.1184 + 0.0007 
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0.12 
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0.13 



Fig. 3 Comparison with other determinations of the strong 
coupling based on data for thrust distribution. The band 
shows the global world average [55]. 



Fig. 2 Scatter plot for the simultaneous fit of cts and ao 
using two different matching schemes. 



ao, the size of the resuhing uncertainty decreases. We 
have assumed that above fij — 2 GeV the physical 
running couphng is weh approximated by its perturba- 
tive component, this allows us to perform a consistency 
check on the running of ao- Since ao is a mean value, 
we consider the integral of the physical cts from 2 GeV 
to 3 GeV 

/ dkas(k^) = 3ao(3GeV) - 2ao(2GeV) = 0.214, 

J 2 GoV 

(44) 

and we compare it to the perturbative result with four 
active flavours 



3GoV 



dkQ^J{k') = 0.218: 



2GcV 



-0.013 
-0.009' 



(45) 



where the values for ao and as{Mz) refer to the find- 
ings in the R matching scheme. We observe that the 
two quantities are in very good agreement within the 
fit uncertainty. This observation indicates that the non- 
perturbative effects on the running of Us can be con- 
sidered negligible above fij = 2 GeV, thereby justifying 
this choice for the non-perturbative matching scale. 

To estimate the impact on as of the finite bottom- 
quark mass corrections we performed a further scan 
without including this correction (setting mf, = 0), us- 
ing therefore only pure massless QCD results. With this 
setting we obtain the following values 



scheme 




ao(2 GeV) 


xVd-o.f. 


R 
log-R 


mil ts+''™22 

U.1110_Q QQJY 

n 1 1 oQ-i-O-0031 

U.iiZO_o 0021 


n c^Qq+0122 

u.ooa_o 

U.OZZ_0.047 


1.56 
1.63 



The values for ao remain almost unchanged while, 
with respect to the nominal scan, the values of as in 
both matching schemes decrease. This can be explained 



by observing that the mass corrections tend to lower the 
thrust distribution in the fit range leading to a higher 
value for a^. The errors on a^. are slightly smaller since 
now we are not accounting for the uncertainty on nib. 



6 Comparison with other otg determinations 

We compare our results to other recent determinations 
of the strong coupling based on thrust. The comparison 
is shown in Fig. 3. The central value is consistent with 
the findings of [31]. On the other hand our theory un- 
certainty is larger. The reduced uncertainty of [31] can 
be explained to some extent by including the sublead- 
ing approximated N'^LL terms, which rely on a Fade 
approximation for the four loop cusp anomalous dimen- 
sion. Such terms carry a scale dependence which leads 
to a reduction of the final error. In our case, we do 
not include the full approximated N'^LL tower, but only 
consider the necessary term for the R matching scheme, 
i.e. Gaia^L, for which an exact analytic expression is 
available [28]. Moreover, a different non-perturbative 
model based on a shape function approach was used 
in [31], and a different technique to perform the renor- 
malon subtraction was implemented. The latter turned 
out to have an important impact on the fit uncertainty. 

In the context of the dispersive model, a global fit 
was performed in [23] to NLL-j-NNLO accuracy in the 
log-R matching scheme, using a larger data set. The 
quoted values for the perturbative coupling and the 
non-perturbative parameter were updated in [56] and 
are as(Mz) = 0.11641°;™^^ and ao(2 GeV) = 0.62 ± 
0.02, respectively. In the log-R matching scheme, we 
note that our resulting fit uncertainty is comparable to 
what was obtained in [56], despite the higher logarith- 
mic accuracy of the resummation in our analysis. This 
is because we consider the simultaneous variation of dif- 
ferent parameters to assess the theory error, leading to 
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a more conservative estimate. The uncertainty in the 
ao value quoted above does not account for the 20% 
variation of the Milan factor. Considering such a sys- 
tematic uncertainty leads to a similar theoretical error 
on qq (c/. erratum to [23]). For completeness, we report 
our fit result obtained with NLL+NNLO perturbative 
accuracy in the log-R scheme: 



scheme 


as{Mz) 


ao(2 GeV) 


xVd-O-f- 


log-R 


U.iiDU_o,o051 


U.OOO_o.o85 1 


1.44 



whose central value is consistent with [23, 56] but 
the uncertainties are larger as expected. 



7 Conclusions 

In the present paper we performed a new fit of the 

strong coupling ag using the event shape observable 
thrust. For the theoretical predictions we use resummed 
NNLL results matched to fixed order NNLO calcula- 
tions. For the hadronisation corrections wc use an an- 
alytical dispersive model, extending the leading renor- 
malon subtraction to 0{al). Furthermore we include 
finite bottom quark mass effects to NLO accuracy. To- 
gether with the perturbative coupling as, a simultane- 
ous fit of the non-perturbative mean effective coupling 
ao was also performed. 

For the as determination we used experimental data 
in the energy range 14 < Q < 209 GeV which have 
been corrected for QED radiation effects from ALEPH, 
L3 and the TASSO experiments. The theoretical un- 
certainties are determined as the envelope of all varia- 
tions obtained by performing a parameter scan in which 
we vary the rcnormalization scale .t^, the resumma- 
tion scale x l , the modified-logarithms parameter p and 
the Milan factor M. . Furthermore we consider bottom- 
quark mass corrections for three different 6-quark masses 
TOfc = 4.0 GeV, rub = 4.5 GeV and = 5.0 GeV. The 
scans are performed using both the R matching scheme, 
leading to the following values 

a«(Mz) =0.113ltnS2i ao(2GeV)= 0.538t°:J°?, 
and the log-R matching scheme which leads to 
a,(Mz) =0.1137t2Z27, ao(2GeV)= 0.524tO;S^t 

In both cases the quoted central value refers to the de- 
fault setting of the theory parameters. The obtained 
values are compatible with recent determinations of as 
from thrust using effective field theory results [31,32], 
but the uncertainty is larger. Compared to previous de- 
terminations of the strong coupling, where Monte Carlo 
hadronisation corrections were used, the present deter- 
mination leads to smaller values for a.. This can be 



mainly explained by observing that Monte Carlo mod- 
els are currently tuned on experimental data using a 
showered (to leading logarithmic accuracy) leading or- 
der perturbative prediction. The effect of the NNLL 
resummation on the determination of has been dis- 
cussed in Section 6, and it shows that higher order loga- 
rithmic terms have a sizeable impact on its final value. 
This suggests that Monte Carlo hadronisation should 
not be used together with higher order resummed pre- 
dictions since this would misestimate the real hadro- 
nisation corrections. This issue has been already dis- 
cussed in [12,25]. Our final value for is also smaller 
than the current world average [55], whose determina- 
tion involves many results from both lattice computa- 
tions and r-decay, which lead to higher values of ag- 
in view of recent developments in extending the resum- 
mation beyond NLL accuracy for event-shape variables 
other than thrust, it would be interesting to repeat a 
similar analysis for the full set of six event-shape observ- 
ables to obtain a more robust fit of the strong coupling 
constant. 
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Appendix A: Explicit resummation formulae 

In the present section we report the full resummation 

formulae used in the main text. The QCD /3-function is 
defined by the renormalisation group equation for the 
QCD coupling constant 

^^ = -a4^)^__/3o + ^/3i + ...j, (A.l) 
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where the first two coefficients read 



The resummation coefficients read 



/3o =^^A - ^Tpnp, 

Pi =%Pa - ^CaTfUf - ^CfTfUf. (A.2) 



Following these conventions, the pcrturbativc functions 
hi{asL) have the following expressions [17,28], as a 
function of A = ^iM^^ log 



hi{asL) 



Ad 



[(1-2A) ln(l-2A) 



A A 

-2(1-A)ln(l-A)], 



(A.3) 



h^iasL) = -:^[21n(l - A) - ln(l - 2A)] 



00 



+ 2-5— ln(l - A) ^[ln(l - 2A) 



1 



hi\l - 2A) - ln(l - A)(2 + ln(l - A)] 



-2^1n.^-^ 



1-2A 



(A.4) 



2B(2) A 



Pi 1-A Pl (1-A)(1-2A) 
2A(2)7J5 A 



+ 



3A2 



PI (1-A)(1-2A) 
ln(l - A) 



(1-A)(1-2A) 



ln(l - 2A) 
1-2A 



1-A 



a 

2^^7s 



Po 



1-A 



2B(i)/3i A + ln(l-A) 



+ 



1 



Pi 1-A ' ^0 (1-A)(1-2A) 

^lEpi 
PI 

(1 - 2A) ln(l - A)] - 7|;A(3 - 2A) 



-[A + (l- A)ln(l-2A) 



.[-A2 + (l-A)(l-2A)lnA__^] 



A^^)pl 



1 - A 



2(1- A)(1-2A) 



ln(l - 2A) 



x[4A + ln(l-2A)]-^3-^^^^^ 

X [A^ - (1 - 2A) ln(l - A)(2A + ln(l - A))] 



(A.5) 



1(1)=": 



AW 



67 
36 



12 



Ca 



18 



(A.6a) 



25 TT- ^ 
324 - 216'+^^''^ 



2051 



7 



1296 72 



/ 15503 389 o 11 , 11 

TT^ H TT* Cs 

V 2592 864 720 4 ^ 



+ CaCf 



48 ^ 36 



Cl 



3 

'32 



17 

96 



^2~ 



3(3 

4 



(A.6b) 



(A.6c) 

3(3 

2 

(A.6d) 



Additional contributions to the functions gi{asL) (8) 
arise from the function .ff (1, as{Q)) which accounts for 

hard virtual corrections as well as the hard collinear 

J [l,as{^I^Q)j and soft 5(l,as(^)) constant 

functions. The latter functions contribute to the pertur- 
bative logarithmic structure since their strong couplings 

are evaluated at the collinear {\J~^Q) and soft {^^jp-) 
scales, respectively. Their perturbative expansions can 
be found for example in the appendix of [27, 28] . 
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